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Translate this text 

A fluid flow sensor uses Particle 
Image Velocimetry to plot the 
vector flow (20) of fluids such as 
the charge in an engine cylinder 



with data processing to determine 
the real center of turbulent 
rotation (38). 
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Die folgenden Angaben sind den vom Anmelder eingereichten Unterlagen entnommen 

) Verfahren zum Messen von Stromungsverhaltnissen von Fluiden 
) Die Erfindung betrifft ein Verfahren zum Messen von 

Stromungsverhaltnissen von Fluiden, insbesondere einer 

Ladungsmasse in Zylindern einer Brennkraftmaschine, 

wobei die Stromung eine Wirbelstromung, insbesondere 

eine Tumble- und/oder eine Drallstromung ausbildet, und 

aus aufeinanderfolgenden Momentanaufnahmen der 

Stromungsverhaltnisse jeweilige Stromungsvektoren mit 

Stromungsrichtung und Stromungsgeschwindigkeit an 

verschiedenen Orten bestimmt und daraus Kennzahlen K 

der Stromung berechnet werden. Hierbei werden aus den 34 

Momentanaufnahmen jeweilige tatsachliche Mittelpunk- 

te der Wirbelstromung ermittelt und die Kennzahlen aus 

den Stromungsvektoren an den verschiedenen Orten be- 
rechnet, wobei zur Berechnung fur jeden Ort eines jewei- 

ligen Stromungsvektors ein zugehoriger Ortsvektor ver- 

wendet wird, der vom tatsachlichen Mittelpunkt der Wir- 
belstromung als Ursprung ausgeht. 
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Beschreibung 

10001] Die Erfindung belrilll ein Verfahren zum Messen von Slromungsverhallnissen von Huiden, insbesondere einer 
Ladungsmasse in Zylindern einer Brennkraftmaschine, wobei die Stromung eine Wirbelstromung, insbesondere eine 
5 Tumble- und/oder eine Drallstromung ausbildet, und aus aufeinanderfolgenden Moment anaufnahmen der Stromungs- 
verhaltnissc jcwciligc Stromungsvcktorcn mil Stromungsrichtung und Slromungsgcschwindigkcit an verschicdcncn Or- 
ten bestimml und daraus Kennzahlen K der Stromung berechnet werden, gemaB dem Oberbegriff des Anspruchs 1, 
[0002] In verschiedenen technischen Gebieten ist es erstrebenswert, eine Fluidbewegung in einem abgeschlossenen 
Raum zu beobachten und zu analysieren. So ist beispielsweise die Lariungsbewegung in einem Brennraum b/w. Zylinder 
10 sowohl fur Otto- als auch Dicselmotorcn ein Parameter, der die Verbrcnnung und dainit die Wirtschaftlichkcit und Lei- 
slungsfahigkeil wesenllich beeinlluBl. Grundsatzlich wird z wise hen Drall, einer Ladungsbewegung um die Zylinder- 
achse, und Tumble, einer Wirbelstromung um eine Achsc senkrecht zur Zylinderachse, unterschieden. Meist trelen beide 
I .adungsbewegungen gleich/eitig in vcrschiedener Form mil 

10003] Aus "Motortcchnischc Zcitschrift" MTZG1 (2000) 1, S. 40 ff ist cs bekannl, Drall- und Tumblekcnnzahlcn mit 

15 Ililfe von mechanisch arbeitenden Aufnahmeverfahren zu crmitteln. Zum Erfassen beider Arten der Ladungsbewegung 
werden an einem Prufsland fur sialioniire Stromung integrative Mcthoden angewendet. Nachtcilig ist jedoch, daB mecha- 
nische Hilfsmittel, wie ein Fliigelrad oder MomentenmeBgerat, zur Aufnahme von Stromungsrichtung und Slromungs- 
gcschwindigkcit vcrwcndci werden und dicsc in der zu messenden Stromung sclbsl angcordnct werden miisscn, wodurch 
die Stromung in unerwunschler Weise beeinfluBt wird. 

20 [0004] Aus der EP 0 566 120 Bl ist es bekannt, zur Messung der Stromungsverhaltnisse u. a. in Brennriiumen von 
Brennkraftmaschinen das sogen. PIV-Verfahren (Particle Image Velocimetry) anzuwenden. Das stromende Fluid enthalt 
hicrbci Markicrungspartikcl (tracer partikcls) und das Fluid wird mittcls eincs kontinuicrlichen oder gcpulstcn Lasers an- 
gestrahll. Dies macht die Fluidverteilung, Stromungsrichtung und -geschwindigkeil millels der Markierungspartikel 
sichtbar. Entsprechende Bilder werden aufgenommen und mittels einer Bildverarbeitung ausgewertet. Es werden bei- 

25 spielsweise mil 25 Hz zweidimensionale Bilder in Form von Vektorfeldern aufgenommen und anschlieBend anhand ei- 
ner Abfolgc von mehrcren aufcinandcrfolgcndcn Vektorfeldern die tatsacnncncn Stromungsgcschwindigkcitcn crmittclt. 
Bei dem PIV-Verfahren in Verbindung mit Tumble-Messungen wird das dreidimensionale Slromungsfeld im Zylinder 
auf die zweidimensionale Ebene reduziert, welche im groBten Zylinderdurchmesser vorliegt. Fur die Ermittlung der 
Drallzahl wird ein um 90° gedrehtes MeBverfahren angewandl, wobei das Slromungsfeld senkrecht zur Zylindcrmitten- 

30 achsc ausgewertet wird. Hicrbci miBt beispielsweise cine Kamera senkrecht von unlcn in den Brennraum. 

[0005] Bei der anschlieBenden malhematischen Berechnung eines Tumble- bzw. Drallkennwerles wird jedoch bei alien 
vorgenannten Verfahren auf einem Ortsfesten Mittelpunkt der Stromung aufsetzt, welcher aber nur theorelisch gegeben 
und im tatsachlichcn Bctricb nicht an der angenommcnen Stcllc vorhanden ist. Dahcr weichen die mit den bekanntcn 
Methoden enriittelten Werte zum Teil erheblich von den real vorhandenen Werten ab. 

35 [0006] Der vorliegenden Erfindung liegt die Aufgabe zugrunde, ein verbessertes Verfahren der obengenannten Art zur 
Verftigung zu stellen, wobei die obengenannten Nachteile uberwunden werden und eine noch genauere Messung und 
Charaktcrisicrung von Wirbclstromungcn crziclt wird. 

[0007] Diese Aufgabe wird erfindungsgemaB durch ein Verfahren der o. g. Art mit den in Anspruch 1 gekennzeichne- 
ten Merkmalen gelost. Vorteilhafte Ausgestaltungen der Erfindung sind in den abhangigen Anspriichen angegeben. 
40 [0008] Dazu ist es erfindungsgemaB vorgesehen, daB aus den Momenlanaufnahmen jeweilige tatsachliche Mittel- 
punklc der Wirbelstromung crmittclt und die Kcnnzahlcn aus Jen Stromungsvcklorcn an den verschicdcncn Ortcn be- 
rechnet werden, wobei zur Berechnung fur jeden Ort eincs jeweiligen Stroniungsvektors ein /ugehoriger Ortsvektor ver- 
wendet wird, der vom taisachlichen Mittelpunkt der Wirbelstromung als Ursprung ausgeht. 

[0009] Dies hat den Vorteil, daB eine Bewegung des gesamten Stromungsfeldes in der Betrachlungsebenc unkritisch 
45 und ohnc EinfluB auf das MeBcrgcbnis ist. Das crfindungsgcmaBe Verfahren benotigt kcinen festen Mittelpunkt. 

|0010) Um die zu messende Stromung sclbst so wenig wie moglich zu bccinflusscn ist es besondcrs vortcilhafl, die 
Momentanaufnahmen mittels eines PIV-Verfahrens (Particle Image Velocimetry) herzustellen. 
[0011] Eine Kennzahl der Wirbelstromung wird gemafi folgender Formel berechnet: 
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L: Langc dcs Zylindcrs, 

D: Durchmesser des Zylinders, 

r( Abstand zwischen MeBpunkt i und realcr Achsc der Wirbelstromung, 

7;: Ortsveklor des MeBpunkles i ausgehend von der realen Achse O der Wirbelslromung, 

(0„ p : Winkclgcschwindigkcit der gcsamtcn McGfia'chc, 

(Omoi: Winkelgeschwindigkeil des Motors, 

Q: DurchfluBwert, 

y„ Km : mitllcrc Gcschwindigkcit, folgcnd aus dcm DurchfluBwert Q und Zylindcrdurchmcsscr D, 

v,: Geschwindigkeitsvektor am MeBpunkt i, 

n: Anzahl von MeBpunkten in der MeBflache (Vektorfeld), 

N: Anzahl von gemessenen MeBflachen (Vektorfeldem), 

K: Kcnnzahl der Wirbelstromung. 

[0012) Als Kennzahl K wird eine Tumblezahl T, welche eine Wirbelstromung in einem Zylinder einer Brennkraftma- 
schine mit Drehachse senkrecht zur Zylinderachse beschreibt, und/oder eine Drallzahl DZ, welche eine Wirbelstromung 
in einem Zylinder einer Brennkraftmaschine mit Drehachse parallel zur Zylinderachse beschreibt, berechnet. 
[0013] Durch das crfindungsgcmaBc Vcrfahrcn ist die Bcrcchnung einer zusatzlichcn charakteristischen GroBc bzw. 
Kcnnzahl moglieh. namlich die Bcstimmung einer Stabilitatszahl S, welche die Bcwegung dcs gcsamtcn Stromungsfcl- 
des beschreibt und gemaB folgender Formel berechnet wird: 



wobei 
und 

und 

r„=±|>, 

r ; : Abstand zwischen realer Wirbelachse i und idealer Wirbelachse O, 65 

r mcan : milllerer Abstand zwischen realer Wirbelachse i und idealer Wirbelachse O, 

r,i cv : Abweichung Abstand zwischen realer Wirbelachse i und idealer Wirbelachse O, 

r max : maximale Abweichung des Abstands zwischen realer Wirbelachse i und idealer Wirbelachse O, 
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D: Durchmesscr des Zylindcrs 

N: Anzahl von gcmcsscnen McBflachcn (Vektorfcldcrn) 
S: Slabilitatszahl. 

|0014] Weitere Merkmale, Vortcile und vorteilhafte Ausgestaltungen der Erfindung ergeben sich aus den abhangigen 
5 Anspriichen, sowie aus der nachslehenden Beschreibung der Erfindung anhand der beigefiigten Zeichnungen. Diese zei- 
gcnin 

|0015J Fig, 1 eine sehematische Darstellung eines slalioniiren Prufslandes fur einen Zylinderkopf, 

[0016] Fig. 2 eine graphische Veranschaulichung des MeBprinzips des PIV-Verfahrens (Particle Image Velocimetry), 

[0017] Fig. 3 eine beispielhafte Momentaufnahme eines MeBfeld zur Berechnung einer I\imblez.ahl, 

10 [0018] Fig. 4 eine beispielhafte Momentaufnahme eines MeBfeld zur Berechnung einer Drallzahl und 
|0019] Fig. 5 eine beispielhafte Momentaufnahme eines MeBfeld zur Berechnung einer Slabililalszahl. 
[0020] Nachfolgend wird die Erfindung anhand der Ermittlung von Kennzahlen einer Ladungsmasse in einem Kolben 
einer Brennkraftmaschine naher beschrieben. Dies isi jedoch lediglich beispielhaft zu verstehen. Das erfindungsgemaBe 
Verf'ahrcn ist zur Analyse jcglichcr Stromungcn gceignct, die cine Wirbclstrbmung ausbildcn. 

15 [0021] Die Fig. 1 und 2 veranschaulichen das PIV-Verfahren (Panicle Image Velocimetry). Wie insbesondere aus Fig. 
1 ersichtlich, wird ein Zylinderkopf 10 mit EinlaBkanalen 12 auf einem Priifstand montiert, der einen Zylinder 14 mit ei- 
nem feststehenden flachbodigen Kolben (nicht dargestellt) ausbildet. Bei einem konstanten Druckunterschied zwischen 
dem Eintritt der EinlaBkanalc 12 und dem Zylinder 14 wird cine Kanalstromung in den Zylinder 14 cingclcitct. Mittcls 
einer Kamera 16 wird mit einem PIV-Verfahren, welches nachfolgend unter Bezugnahme auf Fig. 2 naher erlaulerl wird, 

20 ein zweidimensionales MeBfeld 20 aufgenommen. In einem Computer 18 werden an verschiedenen Orten dieses MeB- 
feldes die entsprechenden Geschwindigkeitsvektoren, d. h. Betrag und Richtung der Geschwindigkeit, analysiert. 
[0022] Wie aus Fig. 2 ersichtlich, wird mittcls des PIV-Verfahrens das sich untcrhalb des Zylindcrkopfcs 10 ausbil- 
dende Strdmungsfeld 26 gemessen. Ilierzu strahlt ein Laser 22 durch eine Lichtschnittoptik 24 auf das Strdmungsfeld 26, 
sodaB sich ein Lichlschnitt 28 ergibt, der das Strdmungsfeld in einer zweidimensionalen Ebene30, d. h. in einem Beob- 

25 achtungsfeld 30, fur die Kamera 16 sichtbar macht. Die Kamera 16 nimmt dabei das Beobachtungsfeld 30 durch ein Ob- 
jcktiv 32 auf. Fur mehrcrc Zcitcn t wird auf diese Wcisc ein zweidimensionales MeBfeld 20 aufgenommen. Entsprc- 
chende Tumble- und Drallzahlen werden in nachfolgend noch naher beschriebener Weise durch numerische Integration 
der fur verschiedene Zeiten t gemessenen zweidiiucnsionalc MeBteld 20 aus den Geschwindigkeilsvekloren und deren 
I, age berechnet. 

J0 [0023] Fig. 3 veranschaulichl die Berechnung einer Tumblezahl, welehc cine Wirbelstromung mil cine Aehsc senk- 
recht zu einer Langsachse des Zylinders 14 beschreibt. Der Zylinder weisl eine Lange L 34 sowie einen Durchmesser D 
36 auf. Hier/.u wird der Lichlschnitt 28 gemaB Fig. 2 in einer Kbene parallel zur Langsachse des Zylinders 14 ausgerich- 
tct und zwar am Or! des groBtcn Durchnicsscrs des Zylindcrs 14, d. h. die Zylindcrlangsachsc vcrlauft genau in der 
Ebene des Lichtschniltes 28. Ilierdurch wird die Tumblestromung im Zylinder 12 erfaBt, die eine Wirbelstromung urn 

35 eine reale Achse O 38 darstellt, die eine Winkelgeschwindigkeit CO exp (t) aufweisl. I ; ur mehrere Zeiten t werden entspre- 
chende 7.weidimensionaleMeBfelder20 aufgenommen. An verschiedenen Orten i mit Ortsvektor rj(t) undendlicherFlii- 
chc A| wird ein cntsprcchcndcr Gcschwindigkcitsvcktor Vi(t) bestimmt. Aus diescn Geschwindigkeitsvektoren V|(t) wird 
dann eine Winkelgeschwindigkeit (Oj am MeBpunkt i ermittelt. Aus alien Winkelgeschwindigkeiten CO, an den verschie- 
denen MeBpunkten i ergibt sich dann die Winkelgeschwindigkeit C0 cxp (t) der gesamten MeBflache. Das Verhaltnis von 

40 (0 elp (t) zu einer Winkelgeschwindigkeit co mol der Brennkraftmaschine ergibt dann die Tumblezahl. 
[0024] Die Berechnung der Tumblezahl crfolgt nach folgcndcn Formcln: 
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65 L: Lange des Zylinders, 

D: Durchmesser des Zylinders, 

r,: Abstand zwischen MeBpunkl i und realer Achse der Wirbelslromung, 

rj: Qrtsvektor des MeBpunkles i ausgchend von der realen Achse 0 der Wirbelstromung, 
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a clv : Winkelgeschwindigkeit der gesamlen MeBflache, 
<O mo i: Winkelgeschwindigkeit des Motors, 
Q: DurehfluBwert, 

V„, can : mittlere Geschwindigkeit, folgend aus dem DurehfluBwert Q und Zylinderdurchmesser D, 

vj: Geschwindigkeitsvektor am MeBpunkt i, 

n: Anzahl von McBpunktcn in der MeBflache (Vcktorfcld), 

N: Anzahl von gemessenen MeBflachen (Veklorfeldern), 

T: Tumble. 

[0025] Das auBergewcihnliche an dem erfindungsgemaBe Verfahren liegl darin, daB im Gegensatz zu herkommlichen 
Bcrcchnungsvcrfahrcn nicht cine vorbestimmtc, feststchende ideale Drehachse als Ursprung fur die Ortsvcktorcn r j(t) 
verwendet wird, sondem stall dessen derjenige Punkt O 38 in detn jeweiligen zweidimensionalen MeBfeld 20, bei dem 
die reale Drehachse der Tumblestromung das zweidimensionale MeBfeld 20 schneidet. 

[0026] In analogerWeise wird die Drallzahl DZ geniaB der Darsiellung in Fig. 4 ermittelt, so daB zur naheren Hrlaute- 
rung auf die obigc Bcschrcibung der Fig. 3 vcrwicscn wird. Im Unterschicd zur lirmittlung der Tumblczahl T gcmaB Fig. 
3 ist der Lichtschnitt senkrecht zur Langsachse des Zylinders 14 ausgerichtet, so daB eine Wirbclstromung (Drall) des 
Tluides urn cine Achse aufgenommen wird, die parallel zur Zylinderliingsachse vcrlauft. 
[0027] Die Berechnung der Drallzahl erfolgt nach folgenden Formeln: 

nD 2 /4 
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L: Lange des Zylinders, 

D: Durchmesser des Zylinders, 

r;: Abstand zwischen MeBpunkt i und rcalcr Achse der Wirbclstromung, 45 
r {. Ortsvcktor des McBpunktcs i ausgchend von der rcalcn Achse O der Wirbclstromung, 
CO„ p : Winkelgeschwindigkeit der gesamten MeBflache, 
co rao[ : Winkelgeschwindigkeit des Motors, 
Q: DurehfluBwert, 

V mC an: rnitllere Geschwindigkeil, folgend aus dem DurehfluBwert Q und Zylinderdurchmesser D, 50 

v,: Geschwindigkeitsvektor am MeBpunkt i, 

n: Anzahl von MeBpunkten in der MeBflache (Vektorfeld), 

N: Anzahl von gemessenen MeBflachen (Vcktorfcldcrn), 

DZ: Drallzahl. 

[0028] Wiederum wird als Ursprung fur die Ortsvektoren r'i(t) derjenige Punkt O 38 in dem jeweiligen zweidimensio- 55 
nalen MeBfeld 20 gewahlt, bei dem die reale Drehachse der Drallstromung das zweidimensionale MeBfeld 20 schneidet. 
[0029J Zusatzlich wird durch das erfindungsgemaBe Verfahren die Berechnung cincr zusiitzlichcn KcnngroBc moglich, 
niimlich der Slabililatszahl, die die Bewegung des gesamten Slromungsfeldes beschreibt. 
[0030] Die Berechnung der Slabililatszahl S erfolgt nach folgenden Formeln: 
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mem N£t ' 

20 

i{. Abstand zwischcn rcalcr Wirbclachsc i und idcalcr Wirbclachsc O, 
r mean : miulerer Abstand zwischen realer Wirbelachse i und idealer Wirbelachse 0, 
r dev : Abweichung Abstand zwischen realer Wirbelachse i und idealer Wirbelachse O, 
25 r [I1M : maximale Abweichung des Abstands zwischen realer Wirbelachse i und idealer Wirbelachse O, 
D: Durchmcsscr des Zylindcrs 

N: Anzahl von gemessenen Mefiflachen (Vektorfeldern) 
S: Slabilitatszahl. 

[0031] Hierbei entspricht der Ursprung O 40 der idealen Wirbelachse. In der Fig. 5 reprasenlieren die Kreuze die ver- 
30 schicdcncn rcalcn Drchachscn in den zcitlich bcabstandct aufgcnonimcncn zwcidimcnsionalcn McBfcldcrn. Es wird cin 
Mittelwert aller Abweichungen der realen Wirbelaehsen von der idealen Wirbelaehsen gebildel. Das Verhallnis dieses 
Mittelwertes zu einem maximalen Abstand ergibt die Slabilitatszahl S. 

BKZUGSZKICHENLISTE 

35 

lOZylinderkopf 
12 EinlaBkanal 
14 Zylinder 
16 Kamera 
40 18 Computer 

20 zwcidimcnsionalcs McBfeld 
22 Laser 

24 Lichtschnittoptik 

26 Stromungsfeld 
45 28Lichtschnitt 

30 zwcidimcnsionalcn Bbcnc/Bcobachtungsfcld 

32 Objektiv 

34 Liinge L 

36 Durchmcsscr D 
50 38 reale Achse der Tumbleslromung 

40 Ursprung ftir Slabilitatszahl 

DZ Drallzahl 

K Kcnnzahl 

T Tumblezahl 

55 

Patentanspriiche 

1. Verfahren zum Messen von Stromungsverhaltnissen von Fluiden, insbesondere einer Ladungsmasse in Zylin- 
dern einer Brennkraftmaschine, wobei die Stromung eine Wirbelstromung, insbesondere eine Tumble- undVoder 

60 eine Drallstrbmung ausbildet, und aus aufeinanderfolgenden Momentanaufnahmen der Strbmungsverhaltnisse je- 

wcilige Strbmungsvcktorcn mit Strdmungsrichtung und Stromungsgcschwindigkcit an verschicdcncn Orten bc- 
stimmt und daraus Kennzahlen K der Stromung berechnel werden, dadurch gekennzeichnet, daB aus den Momen- 
tanaufnahmen jeweilige tatsachliche Mittelpunkte der Wirbelstromung ermittelt und die Kennzahlen aus den Stro- 
mungsvektoren an den verschiedenen Orten berechnel werden, wobei zur Berechnung fur jeden Ort eines jeweili- 

65 gen Stromungsvektors ein zugehoriger Ortsvektor verwendet wird, der vom tatsachlichen Mittelpunkt der Wirbel- 

stromung als Ursprung ausgeht. 

2. Verfahren nach Anspruch 1, dadurch gekennzeichnet, daB die Momentanaufnahmen mitlels eines PlV-Verfah- 
rens (Particle Image Velocimetry) hergestellt werden. 
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3. Verfahren nach Anspruch 1 oder 2, dadurch gekennzeichnet, daB cine Kcnn/ahl dor Wirbclstroinung genial.! fol- 
gender Formcl bercchnet wird, 



V mean 



L: Liinge des Zylinders, 

D: Durchmesser des Zylinders, 

r,: Absland zwischen MeBpunkt i und realer Achse der Wirbelstromung, 

Ortsvcktor des McBpunklcs i ausgchend von dcr rcalcn Achsc 0 der Wirbclstroinung, 
(0t xp : Winkelgeschwindigkeil der gesainten MeBdache, 
(On,^: Winkelgeschwindigkeil des Motors, 
Q: DurchfluBwcrt, 

V mcall : mittlere (jeschwindigkeil, folgend aus dem DurchfluBwert Q und Zylinderdurehmesser D, 

\\: (Jeschwindigkeitsvektor am MeBpunkt i, 

n: Anzahl von MeBpunklen in der MeBflache (Vektorfeld), 

N: Anzahl von gcmcsscncn McBflachcn (Vcklorfcldcrn), 

K: Kennzahl der Wirbelstromung. 

4. Verfahren nach einem der vorhergehenden Anspriiche, dadurch gekennzeichnet, daB als Kennzahl K eine Tlim- 
blezahl T, welche eine Wirbelstromung in einem Zylinder einer Brennkraftmaschine mit Drehachse senkrecht zur 
Zylindcrachsc beschrcibt, und/odcr cine Drallzahi DZ, wclchc cine Wirbelstromung in einem Zylinder einer Brenn- 
kraftmaschine mil Drehachse parallel zur Zylindcrachsc beschrcibt, bercchnet wird. 

5. Verfahren nach einem der vorhergehenden Anspriiche, dadurch gekennzeichnet, daB zusatzlich eine Stabilitats- 
zahl S berechnet wird gemaB der Formel, 



\: Abstand zwischen realer Wirbelachse i und idealer Wirbelachse O, 
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Tnma : mittlerer Absiand /.wischen realer Wirbelachsc i und idealer Wirbelachse 0, 

r dcv : Abwcichung Abstand zwischcnrcalcr Wirbelachsc i und idcalcr Wirbelachsc O, 

r, ml : maximale Abweichung des Abstands zwischen realer Wirbelachse i und idealer Wirbelachse 0, 

D: Durchmesser des Zylinders 

N: Anzahl von gemessenen MeBflachen (Vektorfeldern) 
S: Slabilitatszahl. 
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